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Abstract 

Electroproduction in deep inelastic scattering at HERA is studied in a model considering a finite 
sum of gluon ladders, associated with a truncation of the BFKL series. The approach contains the 
bare two gluon exchange and both one and two rungs contributions. The model is fitted to the 
data on the inclusive structure function F- 2 ,{x,Q 2 ) in the region x < 0.025 and 0.045 < Q 2 < 1500 
GeV 2 , with a good agreement. Such a description for a large span in Q 2 is obtained through a 
suitable modeling of the remaining non-perturbative background. 


1 Introduction 


The high energy limit of the photon-proton scattering has been one of the main open questions 
concerning perturbative QCD and there is a great theoretical challenge in describing such process. 
The successful renormalization group approach summing the contributions of order (o; s In Q 2 )) n , 
namely the DGLAP evolution equations jl]], which has described systematically the deep inelas¬ 
tic data (gluon driven) starts to present slight deviations as the energies reached in the current 
experiments have increased [Q]. Although its limitation has been theoretically determined ||, the 
DGLAP approach has enough flexibility to describe virtual photon initiated reactions on both low 
and high virtualities at small x J|]. Despite this fact, non-linear effects to the standard DGLAP 
formalism, associated with parton (mainly gluons) saturation and unitarity corrections [§J, are 
already known and their importance to describe the relevant observables at small x and estimate 
the future measurements is not negligible |J. 

On the other hand, the high energy QCD calculation encoded in the leading logarithmic ap¬ 
proximation (LLA) BFKL formalism 0 summing the contributions (a s ln(l/x)) n is a powerful 
technique to perform predictions to the physical processes where a clear hard scale takes place, in 
which the approximation of fixed strong coupling (instead of running) in some external scale of the 
process is considered. This evolution equation provides the x evolution of the gluon density at small 
x. The main difficulty regarding this approach is the resulting total cross section for the BFKL 
Pomeron exchange that violates the unitarity limit, stated by the Froissart limit atot < hi 2 (s) [||] 
(y/s is the center of mass energy), in hadronic collisions. It is a general believe that this bound is 
also required in j*p reactions. Recently, the next-to-leading (NLLA) BFKL calculation program 
has been accomplished showing that the convergence of the series and its stability is far from clear 
at the moment Q. Nevertheless, the unitarity problem has been addressed in Ref. ]K]], where 
the multiple LO BFKL Pomeron exchanges dominate over the next order corrections and then the 
final scattering amplitude is unitarized. Such an approach gives rise to an evolution equation for 
the parton densities matching several statements of the saturation approaches mentioned above. 

Recently, we have proposed a different phenomenological way to calculate the total cross sections 


in the hadronic sector |Ll|], i.e. proton-(anti)proton collisions, and in deep inelastic scattering [12]. 
As a starting point, in Ref. [|h| the reliability to describe the pp(pp) total cross sections through a 
QCD inspired calculation was shown to be successful. The main point is: at finite total energies, 
the LLA and NLLA summation implies that the amplitude is represented by a finite sum on terms, 



where the number of terms increases like Ins, rather than by the solution of the BFKL integral 
equation. So one perform a truncation in the BFKL series and use it for a phenomenological 
description of the non-asymptotic energy data. The interest in to take the firts terms in the 
completete series in the truncation is related to the fact that the energies reached by the present 
accelators are not high enough to accommodate a big number of gluons in the ladder rungs that 
eventually hadronize. Corrobating this hypothesis, for example the coefficient weighting the term 
~ In 3 s turn out to be compatible with zero considering even the Tevatron data Jk|, in contrast 
with the expected from the complete resummation. This finding suggest that energies reached by 
the present accelerators are not yet asymptotic. 

Moreover, in the present theoretical scenario the NLO corrections to the LO BFKL approach 
were accomplished, which lead to a decreasing of the hard Pomeron intercept, issues about the 
running coupling constant and modifications in the correspondent anomalous dimension (for a 
short review, see |J). Definitive answers concerning the NLO solution are now far from clear, 
imposing severe limits on the possible applicability of those results to the experimental situation 
(see, for example ||Idj] ). Some remaining pathologies has been cured in the generalization of the high 
energy factorization formula in next-to-leading approximation by Ciafaloni and collaborators [15]. 
The improving procedure is the identification of the collinearly-enhanced physical contributions as 
the most important agent of the instability of the BFKL hierarchy. Although of this reasonably 
determined theoretical status, a comparison of the underlying expressions with the deep inelastic 
data like in the case of LO BFKL was not performed yet. Indeed, currently the main focus is for 
processes with a predominantly hard scale, for instance 7 * 7 * or forward-jets, instead a two scale 
process like DIS. More specificallly, in the full resummation of the BFKL series even at LLA there is 
not enough flexibility to perform modifications to consider the complete subleading contributions. 
A partial calculation of the subleading correction at the LO BFKL has been performed ]I(|, taking 
into account dominant non-leading effects which come from the requirement that the virtuality of 
the exchanged gluons along the gluon ladder is controlled by their transverse momentum squared, 
restricting the available phase-space of the emitted gluons. Such a procedure can be used to 
perform predictions to observables to be measured. 

Regarding the important question of take into account subleading contributions in a fixed order 
expansion of the BFKL approach, this issue was firstly addressed by Fiore et al. in the hadronic sec¬ 
tor [13]. The coeficients of different powers of ln(s) in the series refer to the dominant contribution, 
at asymptotic energies, for each perturbative order. However, when we perform a phenomenolog¬ 
ical description of data the subleading contributions are absorbed in the normalization constants 
for each ladder in the expansion or in the non-perturbative parameters of the model. This is the 
main factor of reliability lying in the data description of Refs. 10 , |l3| and in our more detailed 
study for pp(pp) reactions in Ref. [11]. 

Motivated by the good result of the truncation in the hadronic case, we performed a more 
detailed calculation for those observables and extended the analysis to the non-forward region. 
Concerning high energy limit, the BFKL approach was a natural choice, since it takes into account a 
n-rung ladder contribution. To address a finite phase space for the gluon emission, disregarding the 
resummation and NLO effects intrinsic to the calculations, we consider that the possible exchanged 
ladders are builded by a little number of gluon rungs. A non-asymptotic cross section is obtained 
considering only the one rung contribution (a tot. ~ In s), which is enough to describe the available 
data 0]- The non zero momentum transfer calculation was performed and the role played by a 
suitable choice of the proton impact factor was determined. Moreover, an useful parametrization 
for the elastic slope B e i(s ) [11], consistent with the Regge phenomenology, was introduced allowing 
to describe with good agreement the differential cross section. The main features coming from that 
calculation are also corroborated by the phenomenological studies in Ref. [17]. 

Focusing on the deep inelastic scattering, the similar motivation for the truncation of the 
perturbative series is the small phase space to allow an infinite gluon cascade in the final state. 
Some authors even advocate that the resummation technique in deep inelastic is not completely 








correct [18]. Indeed, in the available energies there is no room on pseudorapidity to enable a 
completely resummed n-rung ladder and studies have reported a strong convergence of the BFKL 
series considering few orders in the expansion (for example, in heavy vector meson production [18]). 
Furthermore, there is important evidence that the asymptotic solution to the BFKL equation is 
inapropriated in the most of the HERA range [18], and the expansion order by order allows to 
identify the onset of the region where the full BFKL series resummation is required. 

In our previous work on this subject [|l^] , we summed up two terms of the perturbative expan¬ 
sion and obtained the imaginary part of the DIS amplitude, hence the inclusive structure function 
F 2 (x,Q 2 ). Using the most recent HERA data on the F^-logarithmic slope to determine the ad¬ 
justable parameters at the small x region, we performed a broadly description of the structure 
function and its slope in the kinematical range of momentum fraction x < 1CU 2 . The resulting F 2 
and its gluon content turn out having a mild logarithmic growth as x diminishes in contrast with 
the steeper increasing from the usual DGLAP or BFKL dynamics, which we discuss a little more 
in the next section. It was also determined that the non-perturbative background is not negligible 
in all the kinematical region and should be better addressed. A remaining question in those previ¬ 
ous works is the role played by the next order in the perturbative expansion and the modeling of 
the soft background, object of this work. To address the issues above, here one performs a more 
detailed study on the inclusive structure function considering the additional two rung contribution 
for the truncated series. With the introduction of this new term the virtual photon-proton cross 
section reaches the behavior settled by the Froissart bound, i.e. gJJ’ ~ ln 2 (l/x), however keeping 
in mind that it is a result coming from a non-asymptotic series. Moreover, generally speaking 
the BFKL Pomeron is only a perturbative approximation to the true Pomeron (valid at a limited 
kinematical range), of which the exact properties are unknown making it necessary to include some 
contribution of non-perturbative physics. We included the non-perturbative Pomeron in two quite 
different forms, described latter on. As a result, the experimental data on the structure function F 2 
are successfully described in a broader region of (Q 2 ,x) variables and a more refined determination 
of the non-perturbative background is performed. 

The final picture is quite similar to the two-pomeron one, where the hard and soft pomerons 
play an equally important role in data description. 

This paper is organized as follows. In the next section one introduces the main formulae and 
the most important features for the inclusive structure function, including the two rung contribu¬ 
tion. In the Sec. (3), an overall fit to the recent deep inelastic data is performed based on the 
present calculations and a suitable modeling for the non-perturbative background in two forms are 
introduced. As a byproduct, the resulting slopes are also shown. Finally, in the last section we 
draw our conclusions. 


2 Finite sum of gluon ladders in deep inelastic scat¬ 
tering 


Regarding the deep inelastic scattering reaction, the total cross section for the process 7 * p —> X, 
where X states for all possible final states, is obtained from Optical theorem through the imaginary 
part of the elastic 7 * p —► 7 * p amplitude. In the limit of very high energies the BFKL approach j|] 
is the most natural approach to treat such a process and is considered to compute the correspondent 
cross section [ 22 ]. As already discussed in the introduction, the main trouble with that framework 
are the unitarity violation coming from the LLA Pomeron exchange and the known strong NLO 
corrections. In previous calculations we have proposed the truncation of the complete series in 
order to obtain an amplitude described by a finite sum of gluon ladders [12]. Below, one presents 
the main kinematical variables, the calculation of the amplitude up to the second order in the 
perturbative expansion and the correspondent F 2 (x,Q 2 ) expression. 

Here, we are interested in the high energy region W 2 >> Q 2 , where W is the center of mass 







energy of the system virtual photon-proton and Q 2 is the virtuality of the probe photon. Defining 

q2 

the momentum fraction, i.e. Bjorken variable, as x ~ qv+w 2 > the inequality above implies x « 1, 
setting the small x regime. The proton inclusive structure function, written in terms of the cross 
sections for the scattering of transverse or longitudinal polarized photons, reads as 


F 2 (x,Q 2 ) = 


Q 2 


4-7T 2 




a T {x,Q 2 ) + cr L (x,Q 2 ) 


(1) 


In the asymptotic high energy limit, for photons with polarization A, the cross section is given 
by the convolution of a perturbative kernel, which provides the dynamics of the process, with the 
corresponding impact factors of the interacting particles, 


0A (x,Q 2 ) = 


d 2 k x d 2 k 2 


(2tt) 


ki 2 


k 2 2 


*1 (ki.) F{X) k x , k 2 ) $p(k 2 ). 


( 2 ) 


Clarifying the notation, Q is the color factor for the color singlet exchange and ki and k 2 are 
the transverse momenta of the exchanged reggeized gluons in the t-charmel. The $)( (ki) is the 
virtual photon impact factor (with A = T, L) and <h p (k 2 ) is the proton impact factor. The first 
one is well known in perturbation theory at leading order [^] , while the latter is modeled since in 
the proton vertex there is no hard scale to allow pQCD calculations. 

The kernel I ? (x,ki,k 2 ) contains the dynamics of the process and has been systematically 
determined in perturbative QCD [9]. The main properties of the LO kernel are well known J5J and 
the results coming from the NLO calculations indicate that the perturbative Pomeron can acquire 
very significant subleading corrections || . The most important feature of the LLA BFKL Pomeron 
is the leading eigenvalue of the kernel, leading to a steep rise with decreasing x, Fix) ~ f , , 

y \nl/x 

where e = 4a s In 2 ~ 0.5. Therefore, the inclusive structure function will present a similar growth 
at low x. Hence, the resulting amplitude and consequently the total cross section or structure 
function, at first glance, overestimates the growing of the observable currently measured. 

We proposed an alternative phenomenological way to calculate the observables at current en¬ 
ergies, showing that a reliable description of both proton-(anti) proton and virtual photon-proton 
collisions is obtained considering the truncation up to two orders in the perturbative expansion 
fill , |l2| . In the accelerators domain the asymptotic regime is not reached and there is no room in 
rapidity to enable an infinite n-gluon cascade, represented diagrammatically by the BFKL ladder. 
Moreover, a steep convergence of the LO BFKL series in few orders in the perturbative expan¬ 


sion has been already reported [|T8| and phenomenological studies indicate that such a procedure 
is reasonable at least in proton-proton collision [^, [l?]]. In order to perform the calculations, it 
should be taken into account the convolution between the photon and proton impact factors and 
the corresponding gluon ladder exchanges in each order. The Born contribution comes from the 
bare two gluon exchange and in the leading order the amplitude is imaginary at high energies and 
written for t = 0 as 


A B °rn(W,t = 0 ) = ^pV 

n f 


e ) 


d 2 ki 


ki 4 X 


(kO^ki). 


(3) 


The next order contribution is obtained from the graphs considering the one rung gluon ladder 
and has the following expression in LLA 


ion"ung(^ t = 0) = e f l n (W 2 /W$) 


cr ki d 2 k 2 
ki 1 k 2 1 


$X*(k 1 )if(k 1 ,k 2 )$ p (k 2 ). 


The a s is the strong coupling constant, considered fixed since we are in the framework of the 
LO BFKL approach. As a remark, the running of the coupling constant contributes significantly 
















for the NLO BFKL, since it is determined by subleading one-loop corrections, for example the 
self-energy and vertex-correction diagrams |J]. The typical energy of the process is denoted by Ho, 
which scales the logarithm of energy and takes an arbitrary value in LLA. 

The perturbative kernel it (ki, k 2 ) can be calculated order by order in the perturbative expan¬ 
sion and is encoded by the BFKL kernel if one considers the LLA resummation. In the present case, 
t = 0 and it describes the gluon ladder evolution in the LLA of ln(s) as already discussed above. 
The Pomeron is attached to the off-shell incoming photon through the quark loop diagrams, where 


the Reggeized gluons are attached to the same and to different quarks in the loop [12]. Since the 
transverse contribution dominates over the longitudinal one, hereafter $]) is the virtual photon 
impact factor averaged over the transverse polarizations p3j], 


f 1 dT f 1 d P k i 2 (l “ 2rr')(l - 2 ppl) 

lj “ 2 Jo 2vr Jo 2vr - 2 ~' ■ - 2 —' 


kl 2 pp' + Q 2 pTT ] 


(4) 


where p , t are the Sudakov variables associated to the momenta in the photon vertex and the 


notation r ; = (1 — t) and p 1 = (1 — p) is used, Ref. [221. 

A well known fact is that we are unable to compute the proton impact factor <h p (k 2 ) using 
perturbation theory since it is determined by the large-scale nucleon dynamics. However, gauge 
invariance requires that <hp(k 2 —»• 0 ) —> 0 and then the proton impact factor can be modeled as a 
phenomenological input obeying that limit and takes a simple form 


4?p(k2) — A/p 


k 2 + p 2 


(5) 


where M p is the unknown normalization of the proton impact factor and p 2 is a scale which is 
typical of the non-perturbative dynamics. Furthermore, these non-perturbative parameters can 


absorb possible subleading contribution in each order of the perturbative expansion [ 12 ], 

Considering the electroproduction process, summing the two first orders in perturbation theory 
we can write the expression for the inclusive structure function, whose contributions have been 
already discussed (section V of Ref. |2 ~c 
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In 
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F?(Q 2 ,P 2 ) 


( 6 ) 


where the functions F^ orn (Q 2 , p 2 ) and F!p(Q 2 , p 2 ) (I meaning one-rung contribution) given by 



(7) 


F^>(Q 2 ,p 2 ) = - in t fi2 


+ f + 2 «3) 


where xq gives the scale to define the logarithms on energy for the LLA BFKL approach (is 
arbitrary and enters as an additional parameter) and £(3) = 5D r (V r3 ) ~ 1-202 is the Riemann 
(■-function. 

Connecting the present result to the Regge phenomenology, the truncation of the perturbative 
series reproduces the main characteristics coming from the Regge-Dipole Pomeron model. In the 
Dipole Pomeron the structure function is F 2 (x,Q 2 ) ~ R(Q 2 ) ln(l/x), whose behavior corresponds 
to the contribution of a double j-pole to the partial amplitude of 7 *p —> 7 *p, where R(Q 2 ) is the 
Pomeron residue function. Moreover, the Dipole Pomeron trajectory has unit intercept, ajp(0) = 1, 
and has been used in several phenomenological fits to hadronic sector and HERA data |24|] . While 
in the Dipole Pomeron picture the residue is factorized from the energy behavior in the amplitude, 
the Q 2 dependence is calculated order by order in perturbation theory in our case. 








In Ref. [Hj, we choose to determine the parameters of the model (/i, xq and M p ) from a smaller 
data set, meaning the latest HERA measurements on the (^-logarithmic derivative reported by 
the HI [[HJ and ZEUS Collaboration (preliminary). The reasons for that are the consistent and 
precise measurements of the slope and the additional fact that using the Q 2 derivative we avoid 
contributions from the non-perturbative background depending weakly on Q 2 . The obtained values 
were consistent with the naive estimates, namely p? would be in the non-perturbative domain 
(fi « A.qcd) in both data sets and xo has a value consistent with the Regge limit. The fitted 
expression to the slope described the data with good agreement, producing effectively the same 
linear behavior in In Q 2 considered by the HI Collaboration fitting analyzes. Considering the x 
dependence of the slope, a gluon distribution softer than those coming from the usual approaches, 
xG(x, Q 2 ) ~ x - A |,0 is obtained towards small x. The growth of the structure function shown 
large deviations from the steep increasing present at both LO BFKL series and DGLAP approach, 
where F 2 ~ x~ x . The non-perturbative contribution (from the soft dynamics), mainly at low Q 2 
virtualities, was found not negligible. One estimated that those effects imply correction of ~ 20 % 
in the overall normalization. 

In the next section we address the effects of introducing the next order in the perturbative 
expansion and a more refined parametrization for the non-perturbative background. In order to 
do so, we need to calculate the two-rung contribution in the ladders summation. Following the 


calculations 02|, |23|] , it results 
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where the notation is the same one of the expressions above and such a contribution should be added 
to the previous results in order to perform a fit to the HERA data. Our final fitting expression is 
then given by 


F 2 (x,Q 2 ) = F?°™ + Fy [, 


n(I) 


one—rune 


+ F^ 11 ’ [two —rung] + F| oft [Background]. 


(9) 


Some remarks about the results derived from the truncation of the perturbative series are in 
order. An important issue emerging is the role of the subleading contributions at a fixed order 
expansion of the BFKL approach, firstly addressed by Fiore et al. in the hadronic sector [|l^]. The 
coefficients of different powers of ln(s) in the series refer to the dominant contribution, at asymptotic 
energies, for each perturbative order. However, performing a phenomenological description of data 
the subleading contributions are absorbed into the normalization constants for each ladder in the 
expansion, or in the non-perturbative parameters of the model. This can shed light about the 
strenght of those contributions in the description of the observables. For instance, in Ref. it 
was found that already at the Tevatron energy range, the coefficient weighting the term ln 3 (l/x) 
is approximately equal to zero in contrast with the expectations from the perturbative coefficient 
calculated in the complete series. This is the main feature of reliability in the data description of 
Refs. [ 


17] and in our more detailed study for pp (pp) reactions in Ref. [12]. 


Concerning the gluon content of F 2 , in general grounds the truncated series provides a mild 
logarithmic growth instead a steep one from DGLAP or LO BFKL resummations. In order to 
illustrate the result coming from adding rungs in the ladder, that is summing new terms in the 
truncation, we can make use of the effective exponent of F 2 . One can parametrize the increas¬ 
ing of the structure function (and the gluon content at small-x) in the simple form F 2 (x,Q 2 ) — 
xG(x,Q 2 ) ~ f(Q 2 )x~ Xe ff, where we can calculate the exponent as A e ff = dlnF 2 /dln(l/x). 

The Born two gluon approximation to F 2 , i.e. first term in Eq. (6), provides a constant 
contribution on ln(l/.x) and therefore it gives A e // = 0. Of coarse, this result is rule out by the 










experimental measurements. When one introduces the one-rung piece, turn out F 2 = fo(Q 2 ) + 
fi(Q 2 ) ln(l/x), where /o and f\ are the Q 2 -dependent (logarithmic) factors appearing in Eq. (6). 
The resulting effective exponent is thus X e ff = [/o//i+ln(l/x)] _1 . If we consider the approximation 
/ 0//1 1 at fixed Q 2 , which it seems reasonable from inspection of Eqs. (7,8), one estimates that 

the effective intercept takes values between A e ff — 0.2 at x = ICG 2 and decreasing to A e ff ~ 0.1 
for x ~ ICG 5 . This growth is softer than the typical DGLAP result A e ff — 0.3. Concerning the 
addition of the two-rung piece, the result is less simple; however if we consider the limit case where 
f2(Q 2 ) fi(Q 2 )j2(Q 2 ), one obtains A e ff ~ 2/ln(l/x). Then, the effective exponent would take 

values 0.4 at x = ICE 2 and saturating at 0.2 around x ~ ICE 5 : values quite similar to the mean 
value from DGLAP expectations. The results estimated above have implications for instance on 
diffraction. The xjp dependence of the diffractive structure function, F<P, depends directly on the 
unintegrated gluon distribution F(xip,kT) = k\ d xG{x ^k\~) I d k\ and then one expects a mild 
logarithmic increasing and a non-trivial (logarithmic) kx dependence. 

Now, having the expression for the inclusive structure function at hand, in the next section we 
compare it with the HERA experimental results, determining the adjustable parameters and the 
range of validity for this model. 


3 Fits to the HERA data: choice of background and 
discussion of the results. 

In order to compare the expression obtained to the structure function T^x, Q 2 ), Eq. i) , with 
the experiment, we choose to use the updated HERA data set starting from the smallest available 
Q 2 for small x region x ~ ICR 2 . The BFKL Pomeron as well as its truncation should be valid 
in a specific kinematical interval, for instance the full series [^] can accommodate data ranging 
from 1 < Q 2 < 150 GeV 2 in a good confidence level. For this purpose we need to include a 
soft background, accounting for the non-perturbative content and providing a smooth transition 
to Q 2 = 0. Notice that an extrapolation to the photoproduction region is still lacking in the 
present analysis, although we have a definite expression for that [[l^]. Thus, here we consider 
the electroproduction case and the low virtualities range will be driven by the non-perturbative 
background, modeled through a soft Pomeron. There exists numerous possibilities from novel 
[|27| , [28| . |29|l to previous |l0| models, which give significant contribution in the data description for 
Q 2 < 10 GeV 2 . In our case, we use the model with the most economical number of parameters 


and for this purpose we have first selected the latest version [31] of the CKMT model 
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Q 2 


A(Q 2 )+1 


Q 2 + a 


where the expression for A(Q 2 ) has the following form: 


AQ") — A 0 1 + 


Q 2 Ai 
A 2 + Q 2 


( 10 ) 


( 11 ) 


where A (Q 2 ) is the Pomeron intercept and the remaining parameters are defined in |3(J. Such a 
model considers a soft Pomeron which is a single pole in the complex angular momentum plane 
having a Q 2 -dependent intercept. Although formally it is not a pure Regge approach, it describes 
the low virtuality region with very good agreement (3(]] with a limited number of adjustable pa¬ 
rameters. The dependence on Q 2 of the structure function comes from the Pomeron residue and 
in general is modeled since there is little theoretical knowledge, namely vertex functions and cou¬ 
plings at the amplitudes. The gauge invariance only requires that it should vanish as Q 2 —> 0. 


Particularly, a model-independent way is more preferable, for instance as performed by [27], where 
the residue function is extracted from data and then fitted with a suitable adjusting function. 






Another possibility is to select a Pomeron for the background which has an intercept equal to 1 
and has the form of a non-perturbative truncated log(^jr) series (soft multipole Pomeron) [ 
with the form 


m, 


F| oft (x,Q 2 ) = Q 2 
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Q 2 + d 


+ b ■ In 



+ c • In 2 



, ( 12 ) 


where this choice has a larger number of parameters than the previous background and we 
naively expect a better accommodation to the data. To proceed, we used two models for the 
background to fit the structure function F 2 . The finite sum of gluon ladders is encoded in the Eqs. 
(jb]) and (||), while the soft Pomeron is given by Eqs. (10) and (12). In both cases we have applied 
the factor (1 — x ) 7 to provide the behavior of F 2 at large x—region and being the same for both 
soft and hard Pomeron, extending the applicability of the models to the x —> 1 region. From the 
dimensional-counting rules these threshold correction factors are given by (1 — x) 2n_1 , where n is 
the spectators number (for the Pomeron it is equal 4). Thus, our considerations are consistent 
with this fact and should provide a good description even at large x. 

For the fitting procedure we consider the data set containing all available HERA data for 
the proton structure function F 2 jlfj, [^, ^0|. Notice that the most recent measurements in HI 
and ZEUS are more accurate (stat. error « 1 %) than the previous ones, providing stringest 
constraints to the parameters. For the fit we have used 496 experimental points for x < 0.025 and 
0.045 < Q 2 < 1500 GeV 2 . We selected the overall normalization factor as a free parameter for the 
hard Pomeron contribution, Eq. (| 6 |), defined as J\f = J2 &f-N p , considering four active flavours. 

In the Figs. (l)-(3) we show the resulting fits considering the two distinct backgrounds. The best 
fit parameters of the model are shown in Tables (1, 2). 


J\f 


x 0 

a s 

A 

a 

A 

Ai 

A 2 

X^ 

0.00312 

1.39 

0.251 

0.2 (fixed) 

0.279 

0.579 

0.108 

1.65 

9.68 

1.14 


Table 1: Parameters of the model with the first background, Eq. (9), obtained from the fit. 


A f 

d 2 

x 0 

OL s 

a 

b 

c 

a 

P 

7 

d 

x x 

0.0191 

0.593 

0.148 

0.242 

0.506 

-0.426 

-0.0495 

0.491 

1.69 

0.727 

0.130 

0.94 


Table 2: Parameters of the model with the second background, Eq. (11), obtained from the fit. 


Before to proceed, we discuss how the results presented here compare with the previous ones 
using only two terms of the perturbative expansion and how the play roled by the background is 
affected by adding a new term in the expansion. Considering only the finite sum of gluon ladders 
(hard Pomeron) one obtains the following: using either the one rung ladder [Eq. (||)] as well 
as the two-rung contribution [Eqs. (j;§),@], the model provides the same y 2 /dof = 1.6 in the 
interval 1.2 < Q 2 < 150 GeV 2 , whereas the fit is degraded for a larger interval of virtualities (i.e. 
y 2 /dof ~ 2.5 for 1.2 < Q 2 < 800 GeV 2 ). These results are in agreement with the analysis |25|| 
(limited to a specific HI data set [^8|), corroborating that the BFKL-like models accommodate 
data at virtualities up to ~ 150 GeV 2 [^]. This kinematical limit corresponds to the region where 
DGLAP resummation start to be important, that is the [cc s ln(Q 2 /Qo)] n terms would dominate 
the a s h\{l/x)] n contributions. The resulting y 2 /dof = 1.6 is still formally large in comparison 
with more phenomenological approaches, however it can be justified if we note the little number of 
parameters considered (3 adjustable constants) and desconsidered the effects of large x. The main 
point is that the addition of a new term does not improve strongly the description of data. This 
fact is supported by the finding that the coefficient weighting the term ln 2 (l/x) in the expansion 
is significantly small in conparison with the ln(l/x) term [jl3|[ . 


































Concerning the role played by the background, the procedure used in the previous work [12] 
was somewhat different. We choose the slope data for fitting the parameters of the model for two 
reasons: (a) they are quite precise and are determined from the updated measurements of 
(b) considering a derivative on Q 2 we are avoiding the contributions from the non-perturbative 
background (depending weakly on Q 2 ), which should be added to the final expression. There, 
we have multiplied the resulting Fi by a factor 1.2, in order to rougly simulate the background 
effects. Of course, this is only a naive estimation and leads to a worse description of data. However, 
the main point is that the background is important in data description over all the kinematical 
range and would be of order ~ 20%. Here, we have modeled carefully the soft background based on 
Regge phenomenology for the Pomeron. From the discussion about the resulting effective exponent 
in the previous section, we naively expect that the background in the one-rung case would have a 
larger effect in mimic the measured A e ff than in the two-rung case. Although we are not able to 
produce a reasonable comparison with the previous work due to the different procedure considered, 
the definitive finding is that the background is active in the full kinematical interval, instead of 
being important only in a specific range, and account for a non-negligible contribution to F^- 

Returning to the present analysis, the final result considering the additional soft background, 
is in agreement with data in a good confidence level for a large span of Q 2 . The x 2 values obtained 
select the logarithm-type as the better background. The negative value for the c constant is a 
non-favoured choice despite the better x 2 /dof. Instead, if we consider only a ~ log (Q 2 /x) or 
~ log 2 (Q 2 /x) parametrization one gets a smaller number of parameters, however the x 2 /dof value 
becomes similar to the CKMT-type. 

In Fig. (1) are shown the results for the low Q 2 data, expected to be dominated by the soft 
Pomeron background. Both choices for the soft Pomeron seem to describe these data bins quite 
well, however the CKMT-type one provides a steeper increasing with x than the log-type up to 
~ 0.25 GeV 2 , which comes directly from the respective intercept for each background; above that 
virtuality both have the same behavior on Bjorken x. 

In Figs. (2-3), one verifies that the description is independent of the specific choice for the 
soft Pomeron, as expected, since in that kinematical region the finite sum of ladders dominates. 
At higher Q 2 the description deviates following the different backgrounds: at Q 2 ~ 90 GeV 2 and 
increasing as virtualities get larger. However, the deviation is present in a kinematical region 
where no data is measured and does not allow an unambiguous conclusion. In general grounds, as 
it is seen from the figures and tables the model of finite sum of ladders with the logarithmic-type 
background describes better the entire data set and probably it would be even better for a wider 
interval of x and Q 2 , but the CKMT-type contribution is phenomenologically preferred due to the 
smaller set of parameters considered. 




Figure 1: The inclusive structure function at very small Q 2 virtualities. The solid line corresponds to the 
model with the first background [Eq. (10), (11)], while the dash-dotted line to the model with the second 
background [Eq. (12)]. 


Comparing the present analysis with the available Regge phenomenology, we have a description 
similar to the two-Pomeron approach p7) : the hard Pomeron in our case is given by the finite sum 
of gluon ladders up to the two-rung contribution and the ([^-dependence is completely determined 
from the perturbative expansion, truncated at order aq. An extrapolation to the photoproduction 
region is still lacking, which would be obtained once the impact factor of the photon is provided at 
Q 2 = 0 and possible singularities are regulated. The hard Pomeron couples to each quarks flavour 
with the same coupling, i.e. it is flavour blind. In a similar way as [27], we naively expect roughly 
to describe the charm content F£ 5 taking | of the hard-pomeron contribution (the fraction means 
e 2 /(e 2 + e d + e s))- We verify that the backgrounds play a significant role not only at small Q 2 , but 
in the whole interval considered, while the pieces separately give good results just for a relatively 
narrow x-region. Undoubtedly this hybrid approach is good for all x and Q 2 if we use a proper 
behavior of F 2 for large x, i.e. the use of an additional non-singlet term. 

The model is still comparable with other Regge models, for instance considering the very 
nice comparative phenomenological analysis of Ref. [29]. There, the two pomeron model, the 
soft dipole Pomeron, a modified two pomeron model and a generalized logarithmic Pomeron are 
analyzed in a large span of Q 2 (including photoproduction). At small-x, our data description is 
in agreement with them, having a similar \' 2 /dof (slightly larger for the CKMT background) and 
with a similar number of free-parameters (even smaller, for instance using CKMT background). 
The same conclusions are valid concerning the logarithmic fitting of F 2 taken into account in Ref. 
p8|] and in the analysis of Ref. pl|]. The main advantage of the present approach is that the 
dependence on Q 2 of the hard pomeron is completely determined from pQCD, in contrast with a 
model-dependent parametrization of the residue pomeron function in the Regge case. 

Considering the QCD fit using the full BFKL series p^], our results without a background 
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Figure 2: The inclusive structure function at small and medium Q 2 virtualities. Same notation of the previous 
plot. 














presents a slightly larger y 2 /dof in a same fitting region (1.5 < Q 2 < 150 GeV 2 ). However, the 
good data description coming from [ 251 is obtained to the cost of a quite low value of the coupling 
constant (leading a lower intercept for the BFKL pomeron), maybe it mimics NLO corrections to 
the pomeron intercept. In our case, we have fixed the coupling constant in the reasonable value 
a s = 0.2, consistent with the virtualities where the fit is performed, and have noticied that the 
effect of left it free is almost negligible, for instance when we consider the logarithmic background. 

In order to perform a detailed comparison between the results above we calculated the slopes 
of the proton structure function, which are presented in the next section. 


4 The Logarithmic derivatives 


The slopes of the proton structure function give valuable information concerning the behavior of 
the gluon distribution and the effective Pomeron intercept. A considerably broad region of DIS 
kinematical region in the upcoming accelerator, i.e. THERA, will effectively enables to probe the 
saturation phenomenon and other asymptotic properties, which can be more explicit in deriva¬ 
tive quantities directly dependent on the gluon content of the proton. For this study it is most 
convenient to consider the logarithmic slopes, defined as follows 


Bq(x,Q 2 ) 


dF 2 (x,Q 2 ) 
d In Q 2 ’ 


(13) 


and 


B x (x, Q 2 ) 


dlnF 2 (x, Q 2 ) 
<91n(l/x) 


For numerical calculations of the slopes we used the following expressions (see ]3l|1), 


(14) 


Bq(x,Q 2 ) 


Q 1 

2A Q 2 


F 2 (x, Q 2 + AQ 2 ) 


F 2 (x, q 2 — aq 2 )J , 


(15) 


F 2 (x + Ax,Q 2 ) - Fo (x - Ax,Q 2 )] 

- m^) -J' (16) 

The early indirect measurements of the slopes contained the shortcoming of correlated bins in 
the (x, Q 2 ) plane |ki|]. The most recent determinations have a better statistics and the kinematical 
variables are no longer correlated ||T9|. Therefore, instead of using the previous two types of slopes 
|34|], which are determined refering two variables, and where one is being averaged (one of them is 
bounded by a kinematical constraint), we calculated the following four types: 


B x (x,Q 2 ) = - 


2Ax 


Bx{x,Q fix) 1 Bq(x, Q' fi x ), 

B x {xfixi Q )j Bq(x fixiQ )■ 


Below, we discuss the results for each slope. In the Fig. (4) is shown the slope BQ(x,Q 2 ix ) 
in virtualities ranging from 2 < Q 2 < 100 GeV 2 . The calculated slopes are in agreement with 
data, describing well the x-dependence. Moreover, the results using the CKMT-type background 
he slightly above the logarithmic-type at all kinematical interval. At larger Q 2 the difference 
between them is more evident, however in such a region there are no available data to provide a 
discrimination. 

In Fig. (5) are shown the results for the Bq(xfi X , Q 2 ) for Bjorken variable ranging from 8.10" " 5 < 
x < 0.02. In the regions where data exist both backgrounds hold. The CKMT-type background 
select parameters such that the description is similar to the results using a QCD fit (HI) [pT The 
logarithmic-type one presents a typical feature: there is a turn-over (a bump) in virtuality Q 2 ~ 15 
GeV 2 , in a similar way as the Regge dipole Pomeron [21]. Both coincide starting at x = 3.2 10 -3 . 














In Fig. (6) are shown the results for the B x (xfi x ,Q 2 ) for Bjorken variable ranging from 
7.1CV' 5 < x < 0.019. Such a quantity corroborates a phenomenological Pomeron having In¬ 
dependent intercept. For this slope, practically the same result is obtained considering the two 
models for the background. A significative deviation is verified only at larger x ~ 0.012 . 

In Fig. (7) are shown the slope B x (x. Q'j lx ) for virtualities ranging from 1.5 < Q 2 < 150 
GeV 2 . One verifies that the slope is independend of Bjorken x for x < 0.01, in agreement with the 
experimental measurements and consistent with the HI NLO QCD fit pel]. The two backgrounds 
deviate from each one in the transition region ~ 0.01, suggesting that the large x region would 
be described differently taking distinct backgrounds. In a rough extrapolation, the CKMT-type 
seems to be favoured, whereas the logarithnric-type would overestimate the large x slope values. 


5 Conclusions 

When looking in the region of high energy limit, we are faced with the lack of connection between 
the Regge approach and perturbative QCD in describing the asymptotic behavior of the hadronic 
(and photon initiated) cross sections. The main issue is whether perturbation theory may shed 
light on the origin and the nature of such physics, i.e the Pomeron induced reactions. In this 
work we have studied in detail the application of the finite sum of gluon ladders, associated with a 
truncated BFKL series, for the inclusive structure function considering up to the two rung ladder 
contribution. This provides the asymptotic behavior expected from the Froissart bound. The 
truncated series describes a broad region on Q 2 , however with a formally large y 2 /dof due to the 
limited number of adjustable parameters. A large span in Q 2 is obtained if we consider a non- 
perturbative background, modeled here as a soft Pomeron. The resulting picture is very close to 
the two-pomeron model of Donnachie-Landshoff, with the hard Pomeron settled by the finite sum 
of ladders. As a result, the structure function is described with good confidence level for data lying 
at x < 0.025 and 0.045 < Q 2 < 1500 GeV 2 . The description of data with x > 0.025 is expected to 
hold, since we have introduced the threshold factor, (1 — x) 7 . There is no unambiguous sensitivity 
to the specific choice of the background, although the CKMT-type is preferred due to the reduced 
number of additional parameters. The region where the deviation between the backgrounds is more 
important stays in the very low virtualities ( Q 2 < 1 GeV 2 ), with the logarithmic-type providing a 
better x 2 /dof. 

Using the mentioned result, the updated slopes of the proton structure function can be numer¬ 
ically calculated considering the two backgrounds. In general grounds, both choices for the non- 
perturbative piece seem to be in agreement with the available data, with deviations in kinematical 
regions where there are not measurements to clarify the analysis. For example, the CKMT-type 
provides a description closer to the NLO QCD fits, favoured due to the reduced number of ad¬ 
ditional constants. For instance, it selects parameters leading to the growth with Q 2 at small x 
for the slope {dF-z/d In Q 2 ) x , and the flat behavior on small Bjorken x at fixed Q 2 for the slope 
( dlnF 2 /dlnx)Q 2 , with similar features as in the QCD fits. 
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Figure 3: The inclusive structure function at large Q 2 virtualities. Same notation of the previous plots. 
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Figure 4: The result for the Bq slope plotted as a function of x for fixed Q 2 virtualities compared with 
the latest HI data [19]. The solid line corresponds to model with the first background [Eq. (9)], while the 
dash-dotted line to the model with the second background [Eq. (11)]. 
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Figure 5: The result for the Bq slope plotted as a function of Q 2 (in GeV 2 ) for fixed x compared with the 
latest HI data (l9| . Same notation of the previous plots. 
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Figure 6: The result for the partial derivative B x slope plotted as a function of Q 2 for fixed x compared with 
the latest HI data [^]. Same notation of the previous plots. 
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Figure 7: The result for the B x slope plotted as a function of x for fixed Q 2 virtualities compared with the 
latest HI data [20|. Same notation of the previous plots. 






























































